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I review recent progress in understanding non-perturbative aspects of string theory, quantum gravity and 
non-commutative geometry using lattice methods. 



1. Introduction 

Lattice methods have been very successful deal- 
ing with string theory, quantum gravity and non- 
commutative field theory in the following sense: 
whenever the bosonic theory is known to exist, 
the lattice formulation has 

(1) provided a rigorous definition of the theory, 

(2) a scaling limit where the continuum theory 
can be defined and which can be analyzed by nu- 
merical and/or analytical methods, 

(3) often been superior to a continuum approach 
even from an analytical point of view. 

Examples are 2d quantum gravity and non- 
critical strings with a central charge c < 1 , 
non-commutative Yang-Mills theory (which was 
in fact discovered on the lattice as early as 1983 
[lj§]), and 3d quantum gravity formulated as a 
Turaev-Viro state-sum. There are also exam- 
ples where lattice attempts to formulate a non- 
perturbative theory have failed. But with the ad- 
vantage of hindsight we now understand that the 
message from the lattice was correct even in these 
cases. One such example is the bosonic string 
theory in dimensions d > 1 ||-[|] . Probably there 
exists no non-tachyonic bosonic string theory in 
dimensions d > 1. 

The examples mentioned are non-trivial in 
the sense that they deal with diffcomorphism- 
invariant theories and one would not expect the 
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lattice to be the best way to regularize such theo- 
ries. In fact, the situation seems to be the oppo- 
site: the lattice framework introduced by the use 
of "dynamical triangulations" seems to work very 
well for the above-mentioned bosonic theories. 

However, there are two fundamental areas 
where the lattice approach cannot yet claim suc- 
cess: four-dimensional quantum gravity and su- 
perstring theory. In the first case we do not know 
if there exists a consistent non-perturbative the- 
ory of quantum gravity, formulated entirely in 
terms of the gravitational fields, or if the theory 
has to be embedded in a larger theory like super- 
string theory. However, if a purely bosonic the- 
ory of four-dimensional quantum gravity exists, 
it follows from the remarks above that the lattice 
framework should be ideal for a non-perturbative 
definition. I will discuss a new lattice approach, 
called Lorentzian simplical quantum gravity in 
Sec. f|. In the case of superstring theories, lattice 
formulations have been stalled by the inability to 
implement supersymmetry in the correct way on 
the lattice. While a number of attempts to imple- 
ment worldsheet supersymmetry on random lat- 
tices have failed, it has also been known that the 
chances were better if one used a Green-Schwarz 
formalism. In that case the supersymmetry is 
a space-time supersymmetry rather than a world 
sheet supersymmetry and thus not in direct con- 
flict with an underlying worldsheet lattice. The 
obstacle came in that case from the so-called k- 
symmetry which has not yet been implemented as 
an exact symmetry in the lattice approach, and 
it is unclear if it will automatically be satisfied 
when the lattice spacing is taken to zero. 
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The so-called "second string revolution" which 
pointed to unexpected and sometimes non- 
perturbative connections between the five dif- 
ferent superstring theories and the so-called M- 
theory, led to the quest for a non-perturbative 
definition of these theories. The new matrix 
model approach to M -theory and type IIB su- 
perstrings initiated by Banks, Fischler, Shenker 
and Susskind (the BFSS matrix model and 
by Ishibashi, Kawai, Kitazawa and Tsuchiya (the 
IKKT matrix model, |8|), respectively, were at- 
tempts to provide such a definition. In these dis- 
cretized models supersymmetry was implemented 
in a way which avoided the problem with k- 
symmetry. The "continuuum limits" of the mod- 
els correspond to N, the size of the matrices, go- 
ing to infinity. However, even for finite N the con- 
cept of supersymmetry has a well defined mean- 
ing in these models. In this respect they differ 
from previous attempts to introduce supersym- 
metry and are more promising. 

2. Superstrings on the lattice 

2.1. Theory 

In the so-called Schild gauge the action of a 
type IIB superstring can be written as the follow- 
ing expression (in dimensions 13 = 4,6 and 10): 

s= Jd^Vg(\{x^x^-^r,{x^i,}), (i) 

where {X, Y} denotes the Poisson bracket be- 
tween the variables ^(^1,^2) and ^(^1,^2) and 
r M are suitably defined 7-matrices. The idea be- 
hind the IKKT matrix model is to replace world- 
sheet variables by N x N matrices: 

(2) 

In this way the worldsheet coordinates (£1,^2) 
are mapped into the matrix indices and the 
worldsheet is replaced by the "matrix lattice" 
(«,j). The continuum limit should be obtained 
in the limit N—>oo. As usual when going from a 
classical theory to a quantum theory we replace 
the Poisson brackets by commutators: 

{X,Y} ^ -i[X,Y], (3) 



and in this way the IKKT action for the type IIB 
superstring becomes: 

S = -Tr (i[A^AH 2 + i#yA^>]). (4) 

The corresponding partition function is 

Z w = | dWdX M e- s '*'*^l. (5) 

As described by IKKT, this model has even for 
finite TV an invariance which might be called su- 
persymmetry. The (bosonic part of the) classical 
equations of motion is 

[X fi ,[X )1 ,X v }] = 0, (6) 

and the simplest solution is 

[X*,X?]=0, ty = 0). (7) 

This means that the matrices X?} can be simul- 
taneously diagonalized and the spray of eigenval- 
ues in R w viewed as the points defining classical 
space-time. 

The valleys X v ] = make the existence of 
the integral (^) non-trivial. The convergence of 
the integral for all N in dimensions D = 4, 6, 10 
(these are the dimensions D > 3 where classical 
supersymmetry of the Green-Schwarz superstring 
can be formulated) was established in ||. 

In a non-perturbative closed string theory the 
dimensionality of real space-time should be deter- 
mined dynamically. We would like to see that 
the typical "quantum" matrices X^ dominating 
the matrix integral (|^) are approximately diag- 
onalizable and that the (suitably defined) spray 
of eigenvalues in the large- N limit constitutes a 
four-dimensional manifold: our world. 

In order to make this more quantitative one can 
define the "space-time" uncertainty A as a mea- 
sure of the ten matrices not being simultaneously 
diagonalizable: 

A 2 = l(TrA 2 - max Y^(UX^)i). (8) 

i 

A 2 =0 if the X^'s are simultaneously diagonaliz- 
able. Let us take the U = U max which minimizes 
the RHS of (||). We define the "space-time coor- 
dinates" 

{ x fi)i = (UmaxX fM U rnax \ . (9) 
\ / ii 
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With this definition we can now define the exten- 
sion of space-time R: 



i? 2 = l(TrA 2 } 



12i<j( X i X j) 



N(N — l)/2 



(10) 



as well as the density of the spray of space-time 
points 



When N = 2 we have an SU(2) matrix model and 
it has been proven that ]10| 



p(r) 



r 



-2D+5 



r large. 



(12) 



There exists good arguments in favor of (|l2|) for 
all values of N @|1§]. 

The simplest way to probe the effective di- 
mension of the space-time dynamically generated 
from the distribution of xf 's is to look at the "mo- 
ments of inertia" for such a space-time: 

T ^-M]v-T)E((<-^)(<™^))- ( 13 ) 



This is a 10 x 10 matrix and the principal mo- 
ments of inertia for the distribution of are the 
eigenvalues Ai > • • • > Aio- 

Our four-dimensional world should appear as a 
flat four-dimensional pancake distribution of x^ 
in the large N limit, i.e. the four first eigenval- 
ues of T^ v should separate from the rest such 
that one has a spontaneous breakdown of the 
ten-dimensional Lorentz invariance in the large- iV 
limit. Further, one would like the ratio A/R — > 
in the large N limit, such that it makes sense to 
talk about a classical background. 

The power behavior ( |i2| ) of p(r) is a source 
of ambiguity in defining R and T^ Vl since 
(Tr X 2n ) = oo for sufficiently large n. This is 
in contrast to the situation where p(r) falls off 
exponentially. 

2.2. Numerical results 

Even if the non-perturbative regularization of 
the type IIB supcrstring has resulted in a finite- 
dimensional "lattice" -theory, the lattice being the 
entries of the NxN matrices X^, the theory is 



not well suited for numerical simulations since it 
has fcrmions. One can integrate out the fermions 
from (pf) and obtain 



Zn = 



J ]JdX fl e i Tr [*«> x "] 2 detM(X), (14) 



where M(X) is a (N 2 -l) x (A 2 -l) matrix. For 
generic X it is complex if the dimension D of 
space-time is 6 or 10, and real and positive if 
D = 4. The first question we want to address 
is whether we observe any trace of spontaneous 
breakdown of Lorentz invariance. Analytical ar- 
guments in favor of a dimensional reduction to 
four dimensions have been given ||, using a one- 
loop approximation, but although very encour- 
aging they need to be substantiated. Computer 
simulations, using the same one-loop approxima- 
tion to the action, can be performed if one drops 
the phase of the determinant in (|lj). Let us de- 
fine 



detM v (X) = |detM(X)| e vT(x) . 



(15) 



v = 1 corresponds to the physical situation we 
want to solve, v — corresponds to the situation 
we can simulate on the computer. There exist 
arguments in favor of a dimensional reduction if 
v = oo |l^] . Consequently, if we observe symmetry 
breaking for v — it would be a strong argument 
in favor of symmetry breaking for v = \. In Fig. 
[j] we show the result of the measurement of the 
eigenvalues of in the case of v = 0. For details 
of the computer simulations shown here, see Jl3[ . 
From the results shown in Fig. [l] there is no trace 
of a spontaneous symmetry breakdown. This is 
an indication that the phase of the determinant 
may play a decisive role in a possible symmetry 
breakdown. However, in Jll[] the model was inves- 
tigated in D = 4 where the fcrmionic determinant 
is real and positive and a dimensional reduction 
to one dimension was observed. This result high- 
lights the ambiguity associated with power law 
distributions like (12). Indeed, T^ v as defined in 
dig] ) is divergent if D = 4 and one of the eigen- 
values of diverges. According to jll| this is 
a reflection of the fact that the tail of the distri- 
bution p{r) is caused by aligned one-dimensional 
configurations of 's. It depends on the choice of 
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Figure 1. The 10 eigenvalues of T M „, extracted 
from Monte Carlo simulations of the one-loop 
approximation to ( |l4| ) with M(X) replaced by 
|M(A)|, and plotted as a function of 1/iV. 



observables how sensitive they are to this tail. As 
an example one can define a modified, converging 
T™ w in L> = 4: 



rpnew 



E 



.(16) 



Using the definition (JlJ) one does not see any 
trace of spontaneously symmetry breaking in the 
case of D = 4. This is illustrated in Fig. |[ Also 



terpretation as expectation values of fundamen- 
tal strings in the theory, do not show any sign of 
spontaneous symmetry breaking |fl3|| . 

This leaves the question concerning the sponta- 
neous breaking of Lorentz invariance in the model 
unsettled. More work is needed in order to un- 
derstand whether or not the signal of breaking 
depends on the observables chosen, and if it does, 
how to interpret this. In addition it would of 
course be very desirable if one could perform sim- 
ulations in D — 10, using the full complex deter- 
minant (or more precisely, Pfaffian). It would 
require handling a complex action, a notorious 
problem in computer simulations. Progress in 
dealing with a problem of the kind encountered 
here has been been made recently [[uj . 

2.3. The BFSS matrix model 

The BFSS matrix model was conjectured to 
provide a non-perturbative definition of M- 
theory. Thus it is an eleven-dimensional theory 
where the matrices X(t) depend on time. The 
problem with fermions also exists in this model, 
and it is more difficult to simulate than the IKKT 
model, since it is a matrix chain model, not a sin- 
gle matrix model. A first attempt has been re- 
ported at this conference fl5|j . Finally, extensive 
numerical studies of a mean field approximation 
to the BFSS model have been performed recently 
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Figure 2. The eigenvalues of T™f; w in the case 
D — A, plotted as a function of 1/N 



measurements of Wilson loops, which have the in- 



3. Non-commutative gauge theories on the 
lattice 

Non-commutative Yang-Mills theory was dis- 
covered on the lattice |l|,|| . In the context of the 
type IIB superstrings discussed in the last section 
they occur by expanding around a classical solu- 
tion to (H) different from (]?]) [ pT|JT^ ], namely one 
satisfying 



(17) 



where C M „ are c-numbers. Note that (17) requires 



the matrices Xf} to be infinite-dimensional. Thus 
one has to work directly at the N = oo limit. Non- 
commutative field theories and in particular non- 
commutative gauge theories have recently been 
studied intensively because they appear both in 
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type IIB string theories and in open string the- 
ory (where the gauge theory, living on the D- 
branes where the open strings end, becomes non- 
commutative in the presence of an external so- 
called B^y field). But non-commutative gauge 
theories can be defined and studied independently 
of string theory and were indeed studied before 
they appeared in string theory. 

As already noted when the theories were 
discovered, a natural regularization of non- 
commutative YM theory is provided by the TEK 
lattice model Q by changing from matrices 



to the exponentials Uu 



This has the fol- 



lowing virtues: (1) It respects non-commutative 
gauge symmetry in the same way as ordinary 
gauge symmetry is respected on the lattice. (2) 
It can be formulated for finite N (this is what 
provides the regularization of the theory) and (3) 
it preserves Morita equivalence. 

The TEK model is U(N) gauge theory on a hy- 
percube with twisted boundary conditions, given 
by the partition function 



Ztek{U{N))= /n dc/ M 



(18) 



where Z Vii — q^u^/n g The classical vac- 

uum of the model becomes non-trivial and can 
be expressed in terms of the so-called twist eaters 



U?=T L 



V V 



(19) 



This equation replaces (|T^). Let N = L D / 2 , D 
even. Then the twist eaters can be used to con- 
struct the so-called Weyl map from the matrices 
U 1 ^ to functions U^x) with arguments x on the 
hypercubic lattice ■ This is the inverse map to 
(^]) . By this map matrix multiplication is mapped 
to the star product of functions: 

U* s ^U{x), U i k u ¥ -> w i ( x ) * W 2 (»), (20) 
k 

where the star product * is defined as 



A(x) *B(x) - 



3 2 * e Jy» z »A{x+y)B(x + z), (21) 



and 0fj, u ~ n^v (sec jl9| for details). The partition 
function Ztek(U(N)) is mapped to 



7 {NC) 



(U(l)) 



JJdU^x) e 



>(w) 



(22) 



UUx+v)*Ut(x), (23) 



where e 2 = g 2 N. This relation: 
Z TEK (U(N)) = zf D C) (U(l)), 



(24) 



which states that a U(N) commutative gauge the- 
ory on a hypercube (a 1 D lattice) is equivalent to 
a non-commutative U(l) theory on a L D lattice, 
IP = N, is the simplest example of the exact 
Morita equivalence for lattice gauge theories jl9| . 
Note that the number of degrees of freedom is 
the same: D N 2 for U(N) on the hypercube and 
D L D for U for the non-commutative U(l) the- 
ory on the L D lattice. It illustrates that lattice 
gauge theory provides us with a natural frame- 
work for defining non-commutative field theories 
in a non-perturbative way. 

4. Quantum gravity on the lattice 

4.1. General considerations 

As mentioned in the Introduction the lattice 
formalism of dynamical triangulations seems the 
ideal framework to address higher dimensional 
quantum gravity. Until now we have not been 
successful^. However, it might be that the prob- 
lems encountered so far are more related to Eu- 
clidean gravity than to quantum gravity. 

In two-dimensional quantum gravity the sim- 
plest direct approach to the theory of fluctu- 
ating geometries has been very successful. By 
rotating to Euclidean signature and regulating 
the sum over geometries by introducing the 
reparameterization-invariant lattice cut-off called 
dynamical triangulations, one was able to calcu- 
late generalized Hartle-Hawking wave function- 
als and correlation functions depending on the 



2 see however the contribution by Shinichi Horata at this 
conference for interesting progress. 
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geodesic distance pofl . However, simple general- 
izations to higher dimensions seem not to work. 
While four-dimensional quantum gravity may not 
exist without being embedded in a larger theory, 
this is not true for three-dimensional quantum 
gravity. This led to the suggestion f2l|| , follow- 
ing an old idea by Teitelboim, that Euclidean 
quantum gravity might not be related to "real" 
Lorentzian quantum gravity, and that one should 
only include causal geometries in the sum over 
histories. The geometries which appear in the 
regularized version of such a theory were called 
Lorentzian dynamical triangulations, and each of 
these geometries has a well defined rotation to 
an Euclidean geometry. The opposite is not true: 
there are many Euclidean geometries which can- 
not be rotated to a Lorentzian geometry with 
a global causal structure. However, it implies 
that one can perform the summation over this 
restricted class of geometries in the "Euclidean 
sector" , and rotate back after the summation has 
been done. This is the way we will treat the sum- 
mation over histories in the following. 

In two dimensions one can perform the sum- 
mation over the class of Lorentzian geometries 
explicitly and obtain a theory which differs from 
Euclidean two-dimensional quantum gravity. The 
difference is best illustrated by considering what 
is called the proper-time propagator, where one 
sums over all geometries with the topology S 1 x 
[0, 1], where the two spatial boundaries are sep- 
arated by a proper time T. In the Lorentzian 
theory the spatial slices at a time T" < T are 
characterized by the fact that the spatial topology 
always remains a circle. In two-dimensional Eu- 
clidean quantum gravity similar spatial slices cor- 
responding to constant proper time split up into 
many (in the continuum limit infinitely many) 
disconnected "baby" universes, each having the 
topology S 1 (22). 

In addition three-dimensional quantum grav- 
ity is interesting for the following reason: locally, 
the classical solution is just flat space, or in the 
case of a positive cosmological constant, 3d de 
Sitter space. If one expands around such a clas- 
sical solution in order to quantize the theory one 
finds it is non-renormalizable. Nevertheless we 
know the theory has no dynamical field degrees 



of freedom, but only a finite number of degrees of 
freedom. Thus it can be quantized following dif- 
ferent procedures, e.g. reduced phase space quan- 
tization. However, it remains unclear if anything 
is "wrong" in an approach where one performs 
the summation over fluctuating three-geometries 
and how such an approach deals with the seeming 
"non-renormalizability" of the theory of fluctuat- 
ing geometries. 

Like in two dimensions, also in three dimen- 
sions there will be a drastic difference between 
what we call Euclidean quantum gravity and 
Lorentzian quantum gravity. Loosely speaking, 
because of the restricted class of geometries which 
enters into the sum over histories in the case of 
Lorentzian quantum gravity, the quantum theory 
will be better behaved, and, contrary to the situ- 
ation in a regularized Euclidean quantum gravity 
theory, one can define a continuum limit. Seem- 
ingly, Euclidean quantum gravity theory, as de- 
fined by dynamical triangulations, does not treat 
the conformal factor correctly (see J2^] for a dis- 
cussion.). It is a little surprising since the main 
success of the formalism in two dimensions pre- 
cisely was the correct treatment of the conformal 
mode. The explanation may be related to the fact 
that the conformal mode in higher dimensions is 
also the cause of the unboundedness of the Eu- 
clidean Einstein action. In the Lorentzian theory 
one still has a conformal mode, but the geome- 
tries associated with this mode are less "singular" 
than the geometries one meets in the Euclidean 
theory (see |24|] for details). 

4.2. 3d Lorentzian dynamical triangula- 
tions 

As mentioned above the proper-time propaga- 
tor is a convenient object to study. We choose 
the simplest possible topology of space-time, 5 2 x 
[0, 1], so that the spatial slices of constant proper 
time have the topology of a two-sphere. Each 
spatial slice has an induced two-dimensional Eu- 
clidean geometry. In the formalism of Lorentzian 
dynamical triangulations the space of Euclidean 
2d geometries is approximated by the space of 2d 
dynamical triangulations. This approximation is 
known to work well and in the limit where the lat- 
tice spacing (the length of the lattice links) goes 
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to zero the continuum limit, i.e. quantum Liou- 
ville theory, is recovered. In order to obtain a 
three-dimensional triangulation of space-time we 
have to fill in the space-time between two succes- 
sive spatial slices. This is done as follows: above 
(and below) each triangle at proper time t = na, 
n an integer, we erect a tetrahedron with its tip 
at t + a, a so-called (3,l)-tetrahedron (if the tip 
is at t — a a (l,3)-tetrahedron). Two tetrahedra 
which share a spatial link in the constant-t plane 
might be glued together along a common time- 
like triangle. Remaining free time-like triangles 
with either the spatial link in the constant-time 
slice at t and the tip at t + a, or the spatial link 
at the constant-time slice at t+a and the tip at t, 
are glued together by so-called (2,2)-tetrahedra. 
They have a spatial link both in the constant-i 
slice and the constant t+a slice. (2, 2) -tetrahedra 
can also be glued to each other in all possible 
ways, the only restriction being that if we cut the 
triangulation in a spatial plane between t and t+a, 
the corresponding graph, which consists of trian- 
gles and squares (coming from cutting the (2,2)- 
simplexes) form a graph with spherical topology. 

Summing over all such piecewise linear geome- 
tries with the Boltzmann weight given by the 
Einstein-Hilbert action defines the sum over ge- 
ometries (see [g4j for details). The partition func- 
tion becomes (up to boundary terms) 

Z(k ,k 3 ,T) = J2c k " N "^- k ^ T \ (25) 

T 

where the summation is over the class of triangu- 
lations mentioned, N a (T) denotes the total num- 
ber of vertices and Ns(T) the total number of 
tetrahedra in the triangulation T. fco is inversely 
proportional to the bare inverse gravitational cou- 
pling constant, while is linearly related to the 
cosmological coupling constant. 

4.3. Numerical simulations 

The model (^5|) can be studied by Monte Carlo 
simulations (see |Q for details). There is only 
one phasef]. Let us fix the total three-volume of 

3 In some previous studies we observed a phase transition 
for large fco . This was caused by restrictions on the gluing 
of (2,2)-simplices. We have now dropped these restric- 
tions. 



space-time to be Na, and let us take the total 
proper time T large. One observes the appear- 
ance of an "semiclassical" lump of universe, as 
shown in Fig. |[ As we change the bare gravita- 
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Figure 3. A snapshot of a three-dimensional con- 
figuration of space-time. The vertical axis is the 
spatial volume N 2 , the horizontal axis the proper 
time t. 

tional coupling constant, the time-extent of the 
semiclassical lump decreases. However, all corre- 
lation functions we have studied can be mapped 
into each other by a simple rescaling of time and 
space directions. Further, we observe that a typ- 
ical spatial volume N 2 (t) in the lump and the 
time-extent AT of the lump scale as: 

N 2 ~ N^ 3 , AT - 7V 3 1/3 . (26) 

This justifies the use of the word "semiclassical" 
in the description of the lump. In the computer 
simulations we observe that the center of mass 
of the lump moves around randomly. In addition 
there are fluctuations in the spatial volume. We 
have studied the fluctuations of successive spatial 
volume. The distribution of such spatial volumes 
is very well described by the formula 

P(N 2 (t),N s (t + a)) ~ e ~ c{ko) mm+m*) . (27) 

The constant c(k ) decreases as fco increases 
(i.e. the bare gravitational coupling constant de- 
creases). At the same time one can observe that 
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the total number of (2,2)-simplices decreases, in- 
dicating that the (2,2)-simplices act as glue be- 
tween successive spatial volumes. 

Thus the leading terms in the effective action 
for the spatial volume of the model are given by 

S ef f(N 2 ) =Jdt (^|| + AN 2 (t)). (28) 

This is exactly the classical Lorentzian action for 
the spatial volume in proper-time gauge, thereby 
supporting the semiclassical interpretation of the 
lump. 

4.4. 3d Lorentzian gravity as a matrix 
model 

If we slice our three-dimensional configura- 
tions, not at proper time t, but at proper time 
t + a/2 we will, as mentioned earlier, obtain a 
spherical graph with two types of triangles, com- 
ing from the spatial intersections of (1,3)- and 
(3,l)-tetrahedra, respectively. In addition the 
graph will contain squares coming from the (2,2)- 
tetrahedra. This class of graphs can be described 
by a two-matrix model: 

Z = JdAdBe- NTl ((A 2 +B^(A* + B> h 0ABAB) m 

A and B are N x N Hermitian matrices and the 
spherical graphs are selected in the large- N limit. 
The coupling constants a, (3 can be related to 
the gravitational coupling constants. While this 
matrix model has not been solved, there exists 
another, closely related matrix model where 

A 3 + B 3 -> A 4 + B 4 , (30) 

which has been solved [^5| . This model has a sim- 
ple interpretation in terms of "triangulations" : 
the (3,1)- and (l,3)-tetrahedra are replaced by 
(4,1)- and (l,4)-pyramids in an obvious notation. 
This model, where the building blocks are pyra- 
mids and (2,2)-tetrahedra, is an equally good reg- 
ularization of three-dimensional gravity. Again 
one can work out the relation between the cou- 
pling constants a, (3 and k ,k 3 (see for de- 
tails). The matrix model is defined for suffi- 
ciently small values of a, f3 and has a critical line 
in the (/3, a)-plane where the continuum limit is 
obtained, see Fig. ||. The three dotted curves ap- 
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Figure 4. The phase diagram of 3d Lorentzian 
quantum gravity from matrix models 



proaching the critical line correspond to three dif- 
ferent fixed values of the bare gravitational con- 
stant. The bare cosmological constant k 3 changes 
along the dotted lines to achieve its critical value 
k^(ko) (which depend on ko) at the critical line 
shown in the figure. Fine tuning of k$ to fcf(fco) 
corresponds to a renormalization of the cosmo- 
logical constant and the approach to the infinite 
volume limit of the theory: for k 3 — > k 3 (ko) the 
expectation value of the number of tetrahedra 
or pyramids diverges. On the other hand there 
seems to be no need for a renormalization of kg: 
it simply defines an overall scale for the model. 

If we follow the critical line of the matrix 
model, starting at small values of j3, it cor- 
responds initially to a weak-coupling phase of 
three-dimensional quantum gravity where the 
bare gravitational coupling constant is small (ko 
large). This is the phase we have observed in the 
computer simulation. However, the matrix model 
undergoes a phase transition at a = /3, shown 
by a dot on Fig. |], separating the weak-coupling 
phase from a strong-coupling phase, correspond- 
ing to large values of the gravitational coupling 
constant. In this phase the triangulations cor- 
responding to the spatial slices at t and t + a of 
topology S 2 disintegrate into many (in the con- 
tinuum limit where a — > into infinitely many) 
baby universes each of topology S 2 , connected by 
a web of thin wormholes. This is possible because 
the matrix model admits more general configura- 
tions than were allowed in the computer Simula- 
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tions. The only requirement of the matrix model 
is that the combined graph at t+a/2 is spherical. 
The component coming from the spatial slice at 
t can actually be disconnected (and similar for 
the component coming from the spatial slice at 
t+a). It becomes a dynamical question what hap- 
pens if one allows for such fluctuations. In the 
weakly coupled phase the spatial topology stays 
unchanged as S 2 , but increasing the gravitational 
coupling constant space starts to be torn apart. 
When the gravitational coupling constant is suf- 
ficiently large a phase transition takes place and 
space disintegrates into baby universes, only con- 
nected by thin wormholes (see [Q for details). 
In this way the model provides a concrete real- 
ization of the ideas of Wheeler and Hawking of a 
quantum foam at short distances. 

Defined in a non-perturbative way on the lat- 
tice, three-dimensional quantum gravity reveals 
a rich structure, and, as in two dimensions, the 
model can be analyzed by a fruitful interplay be- 
tween numerical and analytic methods. 

Hopefully, the same will be true in the case 
of four-dimensional Lorentzian quantum gravity 
defined via dynamical triangulations. 
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